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Abstract 

The evolution of many complex systems, including the world wide web, busi- 
ness and citation networks is encoded in the dynamic web describing the 
interactions between the system's constituents. Despite their irreversible and 
non-equilibrium nature these networks follow Bose statistics and can undergo 
Bose-Einstein condensation. Addressing the dynamical properties of these 
non-equilibrium systems within the framework of equilibrium quantum gases 
predicts that the 'first-mover-advantage', 'fit-get-rich' and 'winner-takes-all' 
phenomena observed in competitive systems are thermodynamically distinct 

phases of the underlying evolving networks. 
PACS numbers: 89.75.Hc, 89.75.-k, 5.65+b 
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Competition for links is a common feature of complex systems: on the world wide web 
sites compete for URLs to enhance their visibility |1|], in business world companies compete 
for links to consumers [0] and in the scientific community scientists and publications compete 
for citations, a measure of their impact on the field . A common feature of these systems 
is that the nodes self-organize into a complex network, whose topology and evolution closely 
reflects the dynamics and outcome of this competition The reversibility of the 

microscopic processes that govern the evolution of these systems, such as growth through 
the addition of new nodes, or the impossibility (citation network) and the high cost (business 
networks) to modify an established links give a distinctly far-from equilibrium character to 
the evolution of theses complex systems. Here we show that despite their non-equilibrium 
and irreversible nature, evolving networks can be mapped into an equilibrium Bose gas [0 , 
nodes corresponding to energy levels and links representing particles. This mapping predicts 
that the common epithets used to characterize competitive systems, such as 'winner-takes- 
air, 'fit-get-rich', or 'first-mover-advantage' emerge naturally as thermodynamically and 
topologically distinct phases of the underlying complex evolving network. In particular, we 
predict that such networks can undergo Bose-Einstein condensation, in which a single node 
captures a macroscopic fraction of links. 

Fitness model- Consider a network that grows through the addition of new nodes such 
as the creation of new webpages, the emergence of new companies or the publication of new 
papers. At each time step we add a new node, connecting it with m links to the nodes 
already present in the system. The rate at which nodes acquire links can vary widely as 
supported by measurements on the www |Q , and by empirical evidence in citation [Q] and 
economic networks. To incorporate the different ability of the nodes to compete for links we 
assign a fitness parameter to each node, 1], chosen from a distribution p{i]), accounting for 
the differences in content of webpages, the quality of products and marketing of companies, 
or the importance of the findings reported in a publication. The probability Ilj that a new 
node connects one of its m links to a node i already present in the network depends on the 
number of links, ki, and on the fitness rji of node i, such that 
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Equation (|T|) incorporates in the simplest possible way the fact that new nodes link prefer- 
entially to nodes with higher k [§] (i.e. connecting to more visible websites; favoring more 
established companies; or citing more cited papers), and with larger fitness (i.e. websites 
with better content; companies with better products and sales practice; or papers with novel 
results). Thus fitness [rji) and the number of links (fcj) jointly determine the attractiveness 
and evolution of a node. 

Mapping to a Bose gas- Using a simple parameterization, we assign an energy to each 
node, determined by its fitness rji through the relation 

ei = -^logr/j, (2) 

where /3 is a parameter playing the role of inverse temperature, j3 = 1/T, whose relevance 
to real networks will be discussed later. A link between two nodes i and j with energy 
and ej (e.g. fitnesses rji and rij) corresponds to two non-interacting particles on the energy 
levels Si and ej (Fig. [I|). Adding a new node to the network corresponds to adding a new 
energy level and 2m particles to the system. Of these 2m particles m are deposited on 
the level ej (corresponding to the m outgoing link node i has), while the other m particles 
are distributed between the other energy levels (representing the links pointing to m nodes 
present in the system), the probability that a particle lands on level i being given by (|^. 
Deposited particles are inert, i.e. they are not allowed to jump to other energy levels. 

Each node (energy level) added to the system at time ti with energy is characterized by 
the occupation number ki(t, e^, ti), denoting the number of links (particles) the node (energy 
level) has at time t. The rate at which level acquires new particles is 

dki{ei,t,ti) e ^'^ ki[ei,t,ti) , . 

dt = z, ' 

where Zt is the partition function, defined as 

Z, = j2e-^^^k,{e„t,t,). (4) 



We assume that each node increases its connectivity following a power-law (we demonstrate 
the self-consistency of this assumption later) 

ki{ei,t,ti) = m i^-j , (5) 

where /(e) is the energy dependent dynamic exponent. Since rj is chosen randomly from the 
distribution p{r]), the energy levels are chosen from the distribution g{e) = /3p(e~^^)e~^^. 
We can now determine Zt by averaging over g{e), i.e. 

<Zt> = J deg{e) ^* rftoe"'^'A;(e, t, to) 

= m (1 + 0(r")), (6) 

where 

is the inverse fugacity. Since z is positive for any /5 7^ we introduce the chemical potential, 
/i, as z = e^^, which allow us to write @ and (^ as 

e-/5M = liin (8) 

t-.oo t ^ ' 

Using (1^) we can solve the continuum equation (|^) finding in a self-consistent way solutions 
of form d^), where the dynamic exponent is 

/(e)=e-^(-^). (9) 

Combining (|^ and (^, we find that the chemical potential is the solution of the equation 

I{P,^^) = jdeg{e) ^^^J^_^ = l. (10) 

The system defined above has a number of properties that make it an unlikely candidate for 
an equilibrium Bose gas First, the inertness of the particles is a non-equilibrium feature, 
in contrast with the ability of particles in a quantum gas to jump between energy levels, 
leading to a temperature driven equilibration. Second, both the number of eligible energy 



levels (nodes) and particles populating them (links) increase linearly in time, in contrast 
with the fixed system size employed in quantum systems. Despite these apparent confiicts, 
Eq. (|T0|) indicates that in the thermodynamic limit (t —>■ oo) the fitness model maps into a 
Bose gas. Indeed, since in an ideal gas of volume f = 1 we have 0] 

deg{e)n{e) = 1, (11) 

where n(e) is the occupation number of a level with energy e Equation (|10|) indicates that 
for the inert gas inspired by the fitness model the occupation number follows the familiar 
Bose statistics 

n(e) = --r-^, , (12) 

i.e. the evolving network maps into a Bose gas. Thus the irreversebility and the inertness 
of the network are resolved by the stationarity of the asymptotic distribution, allowing the 
occupation number to follow Bose statistics in the thermodynamic limit t ^ oo. 

Bose-Einstein condensation- The solutions (^, (^ and exist only when there is a 
that satisfies Eq. ([TD|). However, defined in (|TU|) takes its maximum at /i = 0, thus 

when /(/?, 0) < 1 for a given P and g{e), Eq. (|TUp has no solution. The absence of a solution 
is a well-known signature of Bose-Einstein condensation 0, indicating that a finite no{(3) 
fraction of particles condensate on the lowest energy level. Indeed, due to mass conservation 
at time t we have t energy levels populated by 2mt particles, i.e. 

t 

2mt = ^ k{et(^,t, to) = mt + mtI{P, fj,). (13) 
to=i 



When /(/5, 0) < 1, Eq. (|T|) has to be replaced with 

2mt = mt + mtI{(3,n) + no{(3), (14) 

where no(/3) is given by 

!^ = !-/(/?. 0). (15) 
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The occupancy of the lowest energy level corresponds to the number of links the node with 
the largest fitness has. Thus the emergence of a nonzero riodS), a signature of Bose-Einstein 
condensation in quantum gases, represents a 'winner-takes-all' phenomena for networks, the 
fittest node acquiring a finite fraction of the links, independent of the size of the network. 

The mapping to a Bose gas and the possibility of Bose-Einstein condensation in ran- 
dom networks predicts the existence of three distinct phases characterizing the dynamical 
properties of evolving networks: (a) a scale-free phase, (b) a fit-get-rich phase and (c) a 
Bose-Einstein condensate. Next we discuss each of these possible phase separately. 

(a) Scale-free phase- When all nodes have the same fitness, i.e. p{ri) = 5{rj — 1), {g{e) = 
S{e)), the model reduces to the scale- free model introduced to account for the power-law 
connectivity distribution observed in diverse systems, such as the www actor network 



[^,|^, Internet |T0| or citation networks [0]. The model describes a 'first-mover- wins' 
behavior, in which the oldest nodes acquire most links. Indeed, @ predicts that /(e) = 1/2, 
i.e. according to all nodes increase their connectivity as t^/^, the older nodes with smaller 
ti having larger ki. However, the oldest and 'richest' node is not an absolute winner, since 
its share of links, kmaxi^) / (n^t) , decays to zero as t"^/^ in the thermodynamic limit. Thus 
a continuous hierarchy of large nodes coexist, such that the connectivity distribution P{k), 
giving the probability to have a node with k links, follows a power law P{k) ~ k^^ [HHIl- 
Rewiring, aging, and other local processes can modify the scaling exponents or introduce 



exponential cutoffs in P{k) [JTT]-[T^ while leaving the thermodynamic character of the phase 
unchanged. 

(b) Fit-get-rich (PGR) phase- This phase emerges in systems for which nodes have 
different fitnesses and Eq. (p!OD has a solution (i.e. J(/5,/i) = 1). Equation (^ indicates 
that each node increases its connectivity in time, but the dynamic exponent depends on the 
fitness, being larger for nodes with higher fitness. This allows fitter nodes to join the system 
at a later time and surpass the less fit but older nodes by acquiring links at higher rate 0. 
Consequently, this phase predicts a 'fit-get-rich' phenomena, in which, with time, the fitter 
prevails. But, while there is a clear winner, similar to the scale-free phase the fittest node's 



share of all links decreases to zero in the thermodynamic limit. Indeed, since /(e) < 1, 
the relative connectivity of the fittest node decrease as k{emin^'t) / {fnt) ~ t'^(«min)-i_ This 
competition again leads to the emergence of a hierarchy of a few large 'hubs' accompanied 
by many less connected nodes, P{k) following a power-law P{k) ~ k~"' , where 7 can be 
calculated analytically if p{ri) is known. 

(c) Bose-Einstein (BE) condensate - Bose-Einstein condensation appears when ( |T0|) has 
no solution, at which point (H), (P), ([T0|) break down. In the competition for links the node 
with the largest fitness emerges as a clear winner, a finite fraction of particles (no(/?)) landing 
on this energy level. Thus BE condensation predicts a real 'winner-takes-all' phenomena, in 
which the fittest node is not only the largest, but despite the continuous emergence of new 
nodes that compete for links, it always acquires a finite fraction of links (Eq. ([T5|)). 

To demonstrate the existence of a phase transition from the PGR phase to a BE con- 
densate in a network, we consider a formally simple case, assuming that the energy (fitness) 
distribution follows 



Ce\ (16) 



where 6' is a free parameter and the energies are chosen from e G (0,6^02^), normalization 
giving C = {d+i)/€-^ax- For this class of distributions the condition for a Bose condensation 



IS 



0-1-1 rPe-max 

' dx-^<l, (17) 



where emmit) corresponds to the lowest energy (fittest) node present in the system at time 
t. Extending the limits of integration to zero and to infinity, we find the lower bound for 
the critical temperature Tbe = ^/Pbe 

Tbe > e^a^aO + me + 2))-"^'^'\ (18) 

To demonstrate the emergence of Bose-Einstein condensation near the predicted Tbe-, we 
simulated numerically the discrete network model described above, using the energy distri- 
bution (pIBI) . The chemical potential, /x, measured numerically indicates a sharp transition 



from a positive to a negative value (Fig. 2a), corresponding to the predicted phase transi- 
tion between the BE and the FGR phases. The difference between the network dynamics 
in the two phases is illustrated in Fig.2b, where we plot the relative occupation number of 
the most connected node for different temperatures. We find that the ratio kmaxit) / mt is 
independent on time in the BE phase, indicating that the largest node maintains a finite 
fraction of the total number of links even as the network continues to expand, a signature 
of BE condensation. In contrast, for T > Tbe, the most connected node gradually looses its 
share of links, kmax{t)/Tnt decreasing continuously with time. The numerically determined 
phase diagram (Fig. 2b) confirms that the analytical prediction ([T8|) offers a lower bound 
for Tbe- To predict the precise value of Tbe one needs to incorporate the interplay between 
the convergence of emin{t) to zero and the thermodynamics of BE condensation. 

Since real networks have a fixed (T independent) p{ri) fitness distribution, whether they 
are in the BE or FGR phase is independent of T. Indeed, for p(?7) (1 — rj)^ the network 
undergoes a BE condensation for A > A^^; = 1, thus T, introduced during the calculations, 
plays the role of a dummy variable that at the end vanishes from all topologically relevant 
quantities. The existence of Tbe in the numerically studied model (Fig. 2) is rooted in 
our technically simpler choice of defining g{e) to be independent of T, providing a richer 
phase space to demonstrate the existence of a phase transition in a discrete network model. 
However, as the inset in Fig. 2b shows, by changing 6 the phase transition emerges for fixed 
T as well. 

Do the www, business or citation networks represent a Bose-Einstein condensate, or are 
described by the FGR phase without a dominating winner? It is well known that in certain 
markets some companies did gain and maintain an unusually high market share. A much 
publicized case is the continued dominance of Microsoft in the rapidly expanding operating 
systems market, indicating that business networks do develop 'winner-takes-all' phenomena, 
similar to a Bose-Einstein condensate. The situation of the www is more complex: while 
indeed some webpages did capture an unusually high number of links, it is less clear weather 
they maintain a finite fraction of all links as the www grows (as expected if the www is 
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a BE condensate), or will loose market share, leading to the coexistence of a continuous 
hierarchy of large websites, a signature of the PGR phase. While the power-law connectivity 
distribution holding over six order of magnitude support the FGR scale- free phase P,P 



the most recent study involving more that 200 million nodes did identify a run-away node 
of approximately 10^ links W^, an apparently clear and dominating winner. Since p{r]) 
could be explicitly determined from network data, the phase to which various networks 
belong to could be decided in the future. Uncovering the prevailing behavior in these and 
other complex evolving networks (e.g. social networks [|l5l, ecological and transportation 
webs 1^, etc.) is a formidable task that will require careful quantitative studies on real 
networks. The identification of the possible phases offered in this paper might provide the 
appropriate analytical framework for such studies. 
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""IG. 1. Schematic illustration of the mapping between the network model and the Bose gas. 
(a) On the left we have a network of five nodes (continuous circles and lines), each characterized 
by a fitness, rji, chosen randomly from a distribution p{r]). Equation (Q) assigns an energy 
to each rji, generating a system of random energy levels (right). A link from node i to node j 
corresponds to a particle at level and one at ej. The network evolves by adding a new node 
(dashed circle, tjq) at each timestep which connects to m = 2 other nodes (dashed lines), chosen 
randomly following (|^). In the gas this results in the addition of a new energy level (eg, dashed) 
populated by m = 2 particles, and the deposition of m = 2 other particles at energy levels to 
which ?7e is connected to (e2 and £5). The number of energy levels and particles increase linearly 
with time, as t and 2mt, respectively, (b) In the FGR phase we have a continuous connectivity 
distribution, the several highly connected nodes linking the numerous small nodes together. In the 
energy diagram this corresponds to a decreasing occupation level with increasing energy, (c) In 
the Bose-Einstein condensate the fittest node attracts a finite fraction of all links, corresponding 
to a highly populated energy level, and only sparsely populated higher energies. In (b) and (c) 
the energy diagram shows only incoming links, ignoring the default m = 2 particles on each energy 
level corresponding to the outgoing links to simplify the picture. 
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FIG. 2. Numerical evidence for Bose-Einstein condensation in a network model, (a) Choosing 



the energies from the distribution ([iq) with e. 



1, we calculated the chemical potential 



numerically as the network evolved in time, using (^). Due to the wide range over which varies, we 
plot 1^1 on a logarithmic scale. The temperature at which // changes sign corresponds to the sharp 
drop in |/i| on the figure, and identifies the critical temperature Tbe for Bose-Einstein condensation. 
The data are shown for 6 = 1 and for different times (i.e. system sizes) t = 10^ (continuous), 10^ 
(dashed), 10^ (long-dashed), being averaged over 500,100 and 30 runs respectively. The inset 
shows the chemical potential for different values of the exponent 9 in (16), i.e. 9 = 0.5,1.0,2.0, 
indicating the 9 dependence of Tbe- (b) Fraction of the total number of links connected to the 
most connected ('winner') node, kmax / {in^t) ■, plotted as a function of T, shown for m = 2 and 9 = 1. 
The three curves, corresponding to different system sizes recorded at t = 10^, 10^, 10^ indicate the 
difference between the two phases: in the BE phase (left) the fittest node maintains a finite fraction 
of links even as the system expands, while in the FGR phase (right) the fraction of links connected 
to the most connected node decreases with time. The inset shows the {9, Tbe) phase diagram. 



the continuous line corresponding to the lower bound predicted by Eq. (18), while the symbols 
represent the numerically measured Tbe as indicated by the position of the peaks in Inset (a). 
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